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We try to separate the perturbative and non-perturbative contributions to the plaquette of pure SU (3) gauge 
theory. To do this we look at the large-n asymptotic behaviour of the perturbation series in order to estimate 
the contribution of the as-yet uncalculated terms in the series. We find no evidence for the previously reported 
A 2 contribution to the gluon condensate. Attempting to determine the conventional A 4 condensate gives a value 
~ 0.03(2) GeV 4 , in reasonable agreement with sum rule estimates, though with very large uncertainties. 



1. INTRODUCTION 

There is a long history of trying to extract a 
non-perturbative gluon condensate from lattice 
calculations of the plaquette Q . Thanks to the 
work of [pi we have a much longer perturbation 
series for the plaquette than for any other lat- 
tice quantity. On the other hand, the plaquette 
is one of the most ultra-violet dominated quanti- 
ties, which means that any non-perturbative con- 
tribution is likely to be very small compared with 
the perturbative part. 

Conventionally one expects that the non- 
perturbative part of the plaquette is proportional 
to a 4 , and related to the gluon condensate intro- 
duced in H 
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where bo is the first coefficient of the (3- 
function. However one study has reported a non- 
perturbative contribution scaling like a 2 . 

The Monte Carlo quenched plaquette, Pmc, 
can be measured very accurately. The diffi- 
cult part in determining the condensate is find- 
ing the sum of the perturbative series P per t = 
1 — J2 n P n 9 2n - (We use g 2 as our expansion pa- 
rameter, rather than /3 _1 .) We know the first 
few terms in this series from conventional lat- 



tice perturbation theory calculations. Di Renzo 
et al. have managed to calculate many more terms 
in the series by a clever stochastic method [||. 
Since this is the only lattice quantity with a long 
perturbation series, it gives us a unique opportu- 
nity to see how the coefficients behave at large n. 
The calculations, which were done on a 24 4 lat- 
tice, reproduced the three known terms correctly. 

2. EXTRAPOLATING THE SERIES 

To estimate any non-perturbative contribu- 
tions to the plaquette, we need to know the sum 
of the perturbation series to one part in 10 3 or 
better. We quickly see that for interesting cou- 
plings (g 2 w 1) 10 loops is not enough, and that 
we need a good estimate of the contribution from 
higher order terms. 

To extrapolate to higher loops, we need to know 
how the coefficients depend on n. A good way 
of determining this behaviour is to look at r n , 
the ratio of adjacent coefficients. In statistical 
mechanics, r n is often plotted against 1/n rather 
than n. This is because this gives a straight line 
for a series with a power-law behaviour, 
(1 - ux) q = (2) 
r(n - q) 
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Figure 1. The ratio of successive coefficients in 
the perturbation series of the plaquette, plotted 
against 1/n. The solid line is the fit of eq.^j. 



c„/c„_i = u [1 - (1 + q)/n] 
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where the c n are the coefficients in the Taylor 
expansion. 

The plot, Fig. [l], is almost linear, though there 
is a little curvature, which is taken into account 
by making a fit of the form 
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with u = 0.961(9), q = 0.99(7), s 
This series converges for \g 2 \ < \u[ 
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0.44(10). 
At first 

sight it might seem surprising that a series with 
a power-law singularity describes the data well. 
However, it is known that the specific heat of lat- 
tice gauge theories has a sharp peak in the cross- 
over region between strong and weak coupling (for 
an example, see ||). If this feature dominates the 
series, a fit like eq.Q is reasonable. 

To estimate the higher-order terms in the se- 
ries, [Q use a different extrapolation, based on 
renormalon ideas about the asymptotic form of 
the series. Their formula for the coefficients is 
quite complicated, for full detail, see 0. They 
assume that there is a scheme in which renor- 
malon behaviour (factorial growth of the coef- 
ficients) applies, and that the coupling in this 
scheme is related to the usual lattice coupling by 
a transformation of the type 
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Figure 2. The ratio r n compared with the renor- 
malon fit of ( dashed line ) and the fit of eq. (j^) 
(solid line). 

Using the coefficients known at the time (up to 
n = 8) Burgio et al. made a fit for r and r' which 
leads to the dashed line shown in Fig. |^. Although 
this renormalon fit has about the right value near 
n ~ 8, it does not reproduce the n-dependence as 
well as the fit eq.(Q). 

The resulting fit for the p n is shown in Fig|| 
and compared with the renormalon-inspired fit. 
Eq.(^) describes the coefficients from 2 loops on- 
ward very well, while the renormalon fit shows 
much less curvature than the data. The new 
extrapolation predicts larger coefficients for the 
unknown terms in the perturbation series. That 
means that there will be a smaller condensate left 
after subtracting P per t- 
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Figure 3. The coefficients in the plaquette series, 
compared with the fit of |3/ (dashed line) and of 
eq.fy) (solid line). 
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3. THE GLUON CONDENSATE 

Using the 8 coefficients which they had then, 
and estimating the remainder of the perturbative 
series from their renormalon fit, |Q produced val- 
ues for the non-perturbative part of the plaque- 
tte for (3 between 6.0 and 7.0. The surprising re- 
sult was that AP, the difference between the per- 
turbative plaquette and the Monte Carlo result, 
was not proportional to a 4 A 4 as conventionally 
expected, but to a 2 A 2 . 

In Fig. ^ we compare AP using our estimate 
of the perturbative plaquette, and using the esti- 
mate from [|| . We see that we have a far smaller 
condensate left over, and that there is no longer 
much support for the unconventional a 2 slope. 

In Fig. H we see that AP scales roughly as a 4 , as 
expected from eq. ([!]). However the errors (mainly 
from the uncertainty in summing the perturba- 
tion series) are very large. The dashed lines show 
the values of the gluon condensate expected ac- 
cording to various sum-rule estimates J^Jl . Using 
the force scale r to set our length scale gives 
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which is similar to phenomenological estimates 
such as ~ 0.012 GeV 4 @ or 0.024(8) GeV 4 §]. 

4. CONCLUSIONS 

We can see from Figs that the power law 
fit, eq.(H), describes the known perturbative co- 
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Figure 5. AP, compared with phenomenological 
values of the gluon condensate from 



efficients better than the renormalon fit Q. 

If we use eq.([|) to estimate the sum of the per- 
turbative series, we find a number very close to 
the full plaquette as measured in Monte Carlo cal- 
culations. So any non-perturbative contribution 
to the plaquette is very small. 

The non-perturbative contribution to the pla- 
quette does not look proportional to a 2 A 2 , but 
rather appears to be proportional to a 4 A 4 , with 
about the same order of magnitude as the phe- 
nomenological gluon condensate. 

Asymptotically, the renormalon behaviour 
could still be right. The renormalon asymptote 
predicts r n ss nba/2, which does not catch up 
with eq.(||) until n « 25. 

Acknowledgements: This work was sup- 
ported by the DFG and the BMBF. 



AP 

io~ ; 



3urgio et al. 



Figure 4. AP, the difference between the Monte 
Carlo and perturbative plaquette, according to this 
work (points), and according to /Q/ (dashed line). 
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